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ABSTRACT 


The paper develops a passive sub-millimeter atmospheric 
profile and precipitation retrievals algorithm for Microwave 
Humidity and Temperature Sounder (MWHTS) onboard the 
Chinese Feng Yun 3C (FY-3C) satellite. The retrieval algorithm 
employs a number of neural network estimators trained and 
evaluated using the validated global reference physical model 
NCEP/WRF/ARTS, and works for land and seawater with 
latitude between -30 to 30 degree. To simply the calculation 
procedure and save the training time, principle component 
analysis was adapted to filter out the redundancy caused by 
scanning angle and surface effects, as well as system noise. 
NCEP data per 6 hours are downloaded to run the Weather 
Research and Forecast model WRF, and derive the typical 
precipitation data from the whole world. The Atmospheric 
Radiative Transfer Simulator ARTS is feasible for performing 
simulations of atmospheric radiative transfer. 
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1. INTRODUCTION 


profile; precipitation; 


FY-3C satellite is the new generation polar-orbit meteorological 
satellite of China. FY-3C satellite was successfully launched in 
Taiyuan launch site in Shanxi province of north China, on Sept 
23, 2013 [1-3]. Microwave humidity and temperature sounder 
(MWHTS) onboard FY-3C satellite is a four-frequency, 
fifteen-channel millimeter wave radiometer, which includes the 
mainly sounding channels working at 118.75 GHz for 8 
horizontal polarization channels and 183.31GHz for 5 horizontal 
polarization channels. Also, it has 89GHz and 150GHz (vertical 
polarization) in atmospheric transparent window, shows in Table 
1. Channels working at 118.75GHz are being the first applied 
for polar orbiting satellite to improve the temperature and 
humidity retrieval capabilities of MWHTS, as well as cirrus 
information in the high troposphere. 


Table1. Channel characteristics of MWHTS receivers 


3 | 118.7540 | H 100 3—340 2.0 2.0 
4 | 118.7540 | H 165 3—340 1.6 2.0 
5 | 118.7540 | H 200 3—340 1.6 2.0 
6 | 118.75+1 H 200 3—340 1.6 2.0 
7 | 118.7542 | H 200 3—340 1.6 2.0 
8 | 118.7543 | H 1000 3—340 1.0 2.0 
9 | 118.7545 | H 2000 3—340 1.0 2.0 
1 150.0 V 1500 3—340 1.0 1.3 
1 | 183.3141 H 500 3—340 1.0 1.3 
1 | 183.31+1 H 700 3—340 1.0 1.3 
1 | 183.3143 | H 1000 3—340 1.0 1.3 
1 | 183.3144 | H 2000 3—340 1.0 13 
1 | 183.3147 | H 2000 3—340 1.0 1.3 


Calibrati 
N ieee pola Bandy Dynamic NEAT( on 
0. frequency i zat, idth range (K) K) accurac 
(GHz) ion (MHz) yK) 
1 89.0 V 1500 3—340 1.0 1.3 
2 | 118.7540 | H 20 3—340 3.6 2.0 
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Note: In column 3, V/H means quasi V/H polarization 
2. DATA COLLECTION 


NCEP datasets per 6 hours from January, 1, 2010 to Oct, 31, 
2014 are chosen to run the Weather Research and Forecast 
model WRF [4], and derive the typical precipitation data from 
the whole world, especially for the extreme weather, such as 
typhoon, cyclone and rain form. Bring out the incomplete cases 
and latitude values smaller than 0 degree and larger than 30 
degree. Then among the datasets, the cases which satisfy the 
threshold that precipitation rate larger than 1mm/h are valid. The 
Atmospheric Radiative Transfer Simulator ARTS is used for 
performing simulations of atmospheric radiative transfer, and 
put out the brightness temperature values which respect to 
different angles and surface types. 

In a count, 54 typhoon or cyclone cases are extracted for 
precipitation retrievals training. The whole global datasets 
observed by FY-3C MWHTS from January, 1, 2014 to Oct, 31, 
2014 are used to evaluated and analysis the validation of 
retrieval algorithm. 

For FY-3C MWHTS, channel 9 and channel 15 are sensitive to 
precipitation. Combined with window channels distributed at 
89GHz and 150GHz, people can retrieve precipitation rate using 
neural network estimators. 
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Fig.1. Simulation values of brightness temperature at 118 and 
183GHz (2013.10.01.1600) 


3. DATA PROCESSING 


Comparing simulations from NCEP/WRF/ARTS and FY-3C 
MWHTS observed brightness temperatures for 3 days (about 
2283*8*14 brightness values/day) of data with satellite zenith 
angle distributed at 0°, 5’ ,10° 15° 20° , 30° .40 and 50°. Also, 
MWHTS brightness temperatures have been calibrated against 
global storms simulated by ARTS. 

2) Correction for small brightness temperature biases between 
NCEP/WRE/ARTS and MWHTS 

3) Omission of extreme values 

The following values are omitted: 1) any brightness temperature 
for that footprint is less than 50 K or 

greater than 400 K, which is invalid, 2) the surface altitude is 
above 2 km for |latitude| < 600, or above 1.5 km for 60o 
<|latitude| < 700, or above 0.5 km elsewhere, which could be 
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covered by snow and are sensed more strongly and is called 
too-high, or 3) channel | is less than 240 K, which implies that 
the atmosphere is so cold and potentially dry that precipitation is 


unlikely and that even the most opaque channel, i.e., 183+1 
GHz, may sense the surface and yield false detections of 
precipitation, and is called too-cold. 

4) Principal component analysis (PCA) 

The method of principal component analysis is to extract only 
good signals that provide useful information about precipitation 
but are insensitive to most angle and surface effects, and other 
noises. Only the good signals were used for precipitation 
retrieval. 


4. RETRIEVAL ALGORITHM 


To retrieve precipitation rate for FY-3C MWHTS, neural 
networks are employed to realize training and evaluating for sea 
and land, respectively [5]. As section 3 described, 54 typhoon or 
cyclone cases are simulated by the validated global reference 
physical model, NCEP/WRF/ARTS, composed of the U.S 
National Center for Environment Prediction (NCEP) analyses, 
the new generation National Center for Atmospheric 
Research/Penn State Mesoscale Model (WRF) and the 
Atmospheric Radiative Transfer Simulator, ARTS, which is a 
software for performing simulations of atmospheric radiative 
transfer. 


T profiles 
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Fig 2. Schematic of retrieving flow 


5. RETRIEVAL RESULTS AND 
ANALYSIS 


Microwave humidity and temperature sounder has played an 
important role in the process of monitoring the typhoon. "Fitow" 
typhoon generated on Oct 3, 2013 in the Pacific Northwest 
produce and landed on China in the morning on 16th, 
microwave humidity sounder can achieve high quality data, and 
monitor the typhoon eye area clearly and intuitively. Fig. 4 
shows the brightness temperature of MWHTS, which has the 
resolution of 15 kilometers. In Fig. 4, the blue domain denotes 
the lower brightness temperature values in the typhoon area and 


this is the first time for FY-3C in actual observation applications 
in respect to the typhoon eye area. 

Retrieval results and analysis will be provided in the full paper, 
the author are doing the retrieving work and can make a 
conclusion in the near future. 
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Fig. 3. Retrievals based on brightness temperatures form 
MWHTS 
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Fig. 4. Rain detection of global by FY-3C MWHTS on Jan, 20, 
2015 


For precipitation, the primary radiometric signal at frequencies 
around 183 GHz from precipitating scenes results from the 
scattering by ice hydrometeors [6]. This scattering can result in 
significant brightness temperature depressions (several 10’s K) 
relative to non-precipitating surroundings, and is therefore a 
sensitive proxy for the presence of precipitation at the surface. 
According to the global difference distributions of brightness 
temperatures on Jan 20, 2015, (a) channel 10 plus channel 15, (b) 
for channel 10 plus channel 9 , and (c) for channel 9 plus 
channel 15 for FY-3C MWHTS, using neural network method, 
the precipitation and rain detection can be derived, which is 
shown in Fig.4 and Table 2. 


Table 2. rain detection of typhoon domain by FY-3C MWHTS 
between Jan,1 to Oct,31, 2014 


nesta rain Land sea 
Rate(mm/h) rms/K rms/K 

0.1-1 / 0.84 

1-5 / 3.54 

Typhoon 5-10 / 6.78 
domain 10-30 / 17.21 
30-50 / 21.23 

50-65 / 32.4 


6. SUMMARY AND CONCLUSIONS 


MWHTS plays an important role in studying global climate and 
is the main remote sensing instrument for meteorology and 
disaster. It works in all weather and all day providing the 
observation of brightness temperature which can be used to 
retrieve temperature and humidity profiles and precipitation 
rate. 
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